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ABSTRACT. The algorithm of the ensemble pulsar time scale (PTons) based on the optimal Wiener 
filtration method has been proposed. This algorithm allows the separation of the contributions to the 
post-fit pulsar timing residuals of the atomic clock and pulsar itself. Filters were designed with the use 
of the cross-spectra of the timing residuals. The method has been applied to the timing data of six 
millisecond pulsars. Direct comparison with the classical method of the weighted average showed that 
use of the optimal Wiener filters before averaging allows noticeably to improve the fractional instability 
of the ensemble time scale. Application of the proposed method to the most stable millisecond pulsars 
with the fractional instability CTz < 10^^^ may improve the fractional instability of PTcns up to the level 
~ 10-16. 



1. Introduction 

Despite the algorithms of construction of the ensemble time scales are well established and the frac- 
tional instability of the TAI and TT(BIPM) scales currently is at level of a few IQ-^^ (Petit, 2010), 
development of new methods and approaches of the ensemble scales construction, including pulsar ones, 
presents an interest (Rodin, 2008), (Hobbs et al, 2010). In this paper the method of the optimal Wiener 
filtration is developed for construction of the ensemble pulsar time scale as applied to pulsar timing data. 
This method is compared with the classical algorithm based on the computation of the weighted average 
of time scales. From the general reasoning one can expect increasing of the signal estimation accuracy in 
comparison with the weighted average method since an additional information is used in form of the co- 
variance function or spectrum of the estimated signal. The paper (Rodin, 2008) on the basis of computer 
simulation confirms the above statement. 

In this paper PTens is considered primary as an independent basement for calculation of variations of 
the terrestrial time scales rather than an independent realisation of the barycentric time. 

Section 2 of the paper contains basic formulae of the optimal Wiener filtration. Section 3 tells about 
pulsar observations used in this paper. In section 4 the method of optimal filtration is applied to timing 
data of six millisecond pulsars. The fractional instability is computed for time scales averaged with 
the different methods. 



2. Basic formalae 

Let us assume that we have M time series of the length n ^r(ti) = ^ri{t) {i = l,2,...,n, fc = 
1,2,...,M). In our case ^ri are post-fit timing residuals of time of arrivals (TOAs) of pulsar pulses 
observed relative to the same time scale and with the same registering equipment. This paper uses so 
called square root Wiener filter which expressed by the following formula (Terebizh, 1993) 

where S{uj) is spectrum of the signal, ''N{uj) is spectrum of the noise of the kth pulsar. 

In this paper the problem of stochastic signal estimation is solved under condition of lack of apriori 
information, since the covariance matrix and spectrum of the signal are apriori unknown and estimated 
from the data itself in the assumption that the clock variations (estimated signal) and variations of the 
rotational pulsar phase (additive noise) are uncorrelated. Spectrum S'(w) of the signal were calculated as 



an average of all cross-spectra by the formula (fc 7^ I) 

5M = i-|^-i?M'i?*H|, fc,? = l,2,...,M, (2) 

where (•)* denominates complex conjugation, ^R{ijj) — F\^r{t)\ is Fourier transform of the input data. 
The optimal filter is used according to the following formula 

^R{uj) ^ ^H{ujfR{uj), fc = 1, 2, . . . , M. (3) 

The filtered signal is obtained with the inverse Fourier transform ^s{t) = J^~^[^R{uj)]. 
The averaged signal (the ensemble scale) is calculated by the weighted average formula 

M 

fe=l k 

Weight ''w ^ '^rms~^, rms is the root mean square of the fcth time series. 

3. Observations 

The pulsars PSR J0613-0200, J1640+2224, J1643-1224, J1713+0747, J1939+2134 and J2145-0750 
were observed with the fully steerable RT-64 radio telescope of the Kalyazin Radio Astronomy Obser- 
vatory (KRAO) (Potapov et al., 2003) (Ilyasov et al., 2004), (Ilyasov et al., 2005) (fig. 1). The AS-600 
instrumental facility of the Pushchino Radio Astronomy Observatory (Astro Space Center, Lebedev Phys- 
ical Institute) was used for the registration (Oreshko, 2000). The pulsar pulses were accumulated by a 
spectrum analyzer in two circular polarizations, 80 channels in each, with a channel frequency band 
40 kHz. The observing sessions were conducted, on average, once every two weeks. The total signal 
integration time in each session was about 2 hours. 

The topocentric TOAs were determined by fitting the session-summed pulsar pulse profile into a 
reference template with a hight signal to noise ratio. We computed the barycentric TOAs, determined 
the TOA residulas and refined the pulsar timing parameters by minimizing the residuals by the least- 
squares method using the software Tempo (Taylor, Weisberg, 1989). The DD model (Damour, Deruelle, 
1986) was used to refine and compute the pulsar orbital parameteres. The astrometric, spin and orbital 
parameters of all six pulsars can be found in the paper (Rodin, 2011). 

4. Results 

For the purpose of unification all pulsar data were binned and averaged at the interval 30 days. Gaps 
were filled with the linear interpolation of adjacent values. The common part of the data in the interval 
MJD=51000-53490 were used. Since all observations were carried out with the same registering system, 
no data matching was applied. 

Fig. 2a shows the binned pulsar residuals and their weighted average (solid line). Fig. 2b shows the 
same data passed through the optimal filter ([T]) and their weighted average. 

The stability of a time scale is characterised by so-called Allan variance numerically expressed as a 
second-order difference of the clock phase variations. Since timing analysis usually includes determination 
of the pulsar spin parameters up to at least the first derivative of the rotational frequency, it is equivalent 
to excluding the second order derivative from pulsar TOA residuals and therefore there is no sense in the 
Allan variance. For this reason, for calculation of the fractional instability of a pulsar as a clock, another 
statistic (Tz has been proposed (Taylor, 1991). 

Fig. 3 shows the fractional instability of the ensemble pulsar time scale constructed on the basis of 
six millisecond pulsars observed at KRAO. From the fig. 3 one can see that the optimal filter applied to 
data before weighted average improves the fractional instability almost two times along all time interval 
T. On the basis of visual analysis of the fig. 2 one can conclude that the common signal presented in all 
pulsar data and caused by behavior of the registering equipment and by variations of the local frequency 
standard is better determined in data passed through the filter. 

At the present time at least five pulsars display the fractional instability at the level < 10~^^ (Verbiest 
et al, 2009). One can expect that application to them of the optimal filter method will allow to obtain the 
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Figure 1: Post-fit timing residuals of six millisecond pulsars observed at the Kalyazin Radio Astronomy 
Observatory. 




2 3 4 5 
T-1998.5 



40 

20 

-20 
-40 



(b) 



1 2 3 4 5 6 
T-1998.5 



Figure 2: The binned post-fit timing residuals of six millisecond pulsars observed at the Kalyazin Radio 
Astronomy Observatory before (a) and after (b) applying Wiener optimal filter. Solid line indicates 
weighted average. 
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Figure 3: Fractional instability CTz of the ensemble pulsar time scale in dependence on the observation 
interval r. Dashed line indicates weighted average method, solid line indicates optimal filters method. 



fractional instability of the ensemble pulsar time scale at the level of a few units of 10 , i.e. comparable 
with the instability of the best terrestrial frequency standards. 
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